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Abstract. Using a combined quantum mechanical/classical method, we study the dynamics of deposition 
of small Na clusters on Ar(OOl) surface. We work out basic mechanisms by systematic variation of substrate 
activity, impact energy, cluster orientations, cluster sizes, and charges. The soft Ar material is found to 
serve as an extremely efficient shock absorber which provides cluster capture in a broad range of impact 
energies. Reflection is only observed in combination with destruction of the substrate. The kinetic energy 
of the impinging cluster is rapidly transfered at first impact. The distribution of the collision energy over 
the substrate proceeds very fast with velocity of sound. The full thermalization of ionic and atomic energies 
goes at a much slower pace with times of several ps. Charged clusters are found to have a much stronger 
interface interaction and thus get in significantly closer contact with the surface. 

PACS. 36.40.Gk, 36.40.Mr,36.40.Sx, 36.40.Vz, 61.46. Be Atomic and molecular clusters 



1 Introduction 

Clusters on surfaces are an appealing area of research 
which have motivated numerous studies [1]. The litera- 
ture on the topic is flourishing and many conferences have 
been (fully or partly) devoted to these researches, see for 
example the series of recent ISSPIC conferences [IIEIIII 
[5]. It is now possible to make a direct deposition of size 
selected clusters on a substrate [6l[7]. This opens up new 
possibilities for the synthesis of nano-structured surfaces. 
But the deposition process is not necessarily simple and it 
may lead to a significant modification of the cluster, both 
in terms of its electronic structure and of its ionic ge- 
ometry. There is a subtle interplay of the interface energy, 
electronic band structure of the substrate, and surface cor- 
rugation. These various aspects have already been investi- 
gated in great detail, especially from the structural point 
of view, both experimentally piQlflO] and theoretically [TTl 
fT2 l [T3 l fT4 l [T5 l [T6 l lT7 t fT8]. The situation becomes even more 
involved when one considers the deposition dynamics it- 
self. But that makes the case also more interesting. And 
in a such involved situation, one should start to disen- 
tangle the various influences by studying simple systems, 
i.e. simple geometries and material combinations. In this 
contribution, we are presenting first results for the depo- 
sition dynamics of small Na clusters on Ar(OOl) surface. 
To that end, we employ a hierarchical model which treats 
the different subsystems at different levels of refinement, 
depending on their relevance for the whole process. 

While the presence of a substrate makes the exper- 
imental handling of clusters easier, it strongly compli- 
cates the theoretical description because of the huge num- 



ber of degrees of freedom of the surface. Most theoreti- 
cal approaches thus remain limited to molecular dynam- 
ics (MD). The MD provides the cheapest way to access, 
at least in a gross way, the dynamics of the substrate. 
However, its applicability remains restricted to a narrow 
range of physical situations, because a proper description 
of electronic degrees of freedom is missing. It is never- 
theless crucial to try to account for the surface electronic 
degrees of freedom, in particular if non-adiabatic processes 
become involved. 

A first simplification is to use relatively simple clus- 
ter/substrate combinations as, e.g., simple metal cluster 
on an insulator surface. Because of its inert nature the 
insulator surface can be included at a lower level of de- 
scription. This was, e.g., explored for the case of Na clus- 
ters on NaCl in [TSlfTO] where the cluster electrons were 
described fully quantum-mechanically while the substrate 
is frozen and just serves to deliver an effective interface 
potential (itself tuned to ab-initio calculations [14j). How- 
ever rare-gas atoms show sizeable polarization response. 
These effects need to be included. A somewhat better de- 
scription of surface degrees of freedom can be achieved 
by allowing for a minimum of dynamical response of the 
substrate, as recently proposed for the description of Na 
clusters embedded in rare gases pOll^l^ . The electronic 
degrees of freedom of the cluster were treated microscop- 
ically by time-dependent density-functional theory while 
Na ions and substrate atoms were handled with classical 
MD, taking carefully into account the dynamical polar- 
ization of the atoms. This model belongs to the family of 
coupled Quantum-Mechanical with Molecular- Mechanical 
methods (QM/MM) which are often used in other fields 
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as, e.g., bio-chemistry [23 l [24 l l25| or surface physics [26l 

In this paper we apply this hierarchical approach to 
study the deposition of small Na clusters on an Ar sur- 
face. The paper is organized as follows: Section [2] provides 
a brief presentation of the hierarchical approach and the 
construction of a surface. In section [3l we discuss the ef- 
fects of the dynamical treatment of the Ar atoms. In sec- 
tion m we discuss variation of coUisional energy in terms 
of ionic motion and energetic obscrvablcs. In section [5l 
we study the influence of the initial coUisional geometry 
And in section [6l we look at variation of cluster size and 
charge. 

2 Model 

We give here a short summary of the model, for more de- 
tails see appendix [Al and [28] . The Na cluster is described 
in the TDLDA-MD approach, which was well validated 
for linear and non-linear dynamics of free metal clusters 
p9l[30] . The Na electrons are treated by means of density- 
functional theory at the level of the time-dependent local- 
density approximation (TDLDA), and the Na^ ions are 
propagated by means of molecular dynamics (MD). The 
electron- ion interaction is described by soft, local pseudo- 
potentials [31]. Each Ar atom has two classical degrees 
of freedom: ccnter-of-mass and electrical dipolc moment. 
The dipoles allow an explicit treatment of the dynamical 
polarizability of the atoms. The key constituent of the dy- 
namic cluster-atom interaction are polarization potentials 
[32] . A local pseudo-potential is added for the electron- Ar 
short-range repulsion, modeled following [33], and with a 
final slight readjustment to the NaAr dimer as benchmark 
(bond length, binding energy, and optical excitation). The 
Na-Ar Van-der-Waals interaction is computed via the vari- 
ance of dipole operators [20..28u33j . The atom-atom inter- 
actions are described by a standard Lcnnard-Jones poten- 
tial, while the Ar-Na"*" subsystem is treated by means of 
effective potentials from the literature [34]. More details 
on these various components, including the choice of pa- 
rameters can be found in appendix. 

The total energy thus composed is the starting point 
for variation leading to the time-dependent Kohn-Sham 
equations for the cluster electrons coupled with Hamil- 
tonian equations of motion for the classical degrees of 
freedom (Na"*" ions, Ar positions and dipoles). The ini- 
tial condition is obtained from solving the corresponding 
stationary equations and finally boosting the Na cluster 
to the wanted impact energy. 

The numerical solution proceeds with standard meth- 
ods [30] . We use space-grid techniques to solve the (time- 
dependent) Kohn-Sham equations for the cluster electrons. 
The time propagation is based on a time-splitting method, 
and the stationary solution is attained by accelerated gra- 
dient iterations. The electronic mean-field is treated in ax- 
ially averaged approximation |35[I36] which was found to 
provide an acceptable approximation for the present test 
case of the nearly axial Nag cluster impinging on a sur- 
face. In the following, the symmetry axis will be denoted 



by the z axis. It also corresponds to the deposition axis. 
The Na+ ions and the Ar atoms are, of course, treated in 
full three dimensions. 

Several observables can be computed. We employ here 
simple geometric and energetic indicators, the ionic/atomic 
positions and their kinetic energies. We have checked also 
electronic obscrvablcs, as dipolc oscillations and ioniza- 
tion. Both play a minor role for the analysis (e.g., elec- 
tronic emission during deposition is less than 0.02 %). 
However, this does not imply that the dynamical path 
evolves along a simple and geometrically prescribed Born- 
Oppenheimer surface. Cluster electrons and ions couple to 
the many atoms in the surface through short range colli- 
sion and long-range dynamical polarizability which excites 
all degrees of freedom in a state far from equilibrium. 

Our main test case is Nag on Ar surface. Nag consists 
out of five Na ions in a ring plus one 6th ion topping the 
ring. The top ion sits on the symmetry axis. The Ar(OOl) 
surface is modeled as six layers of 8x8 squares containing 
together 384 Ar atoms. The squares are copied periodi- 
cally in both horizontal directions to simulate an infinite 
surface. To stabilize the underlying (supposedly infinite) 
crystal structure, the atomic positions in the lower two 
layers are frozen at the bulk positions. Some comparisons 
have been performed with the larger substrate configura- 
tion Ar5i2, composed by 6 active layers and 2 frozen ones. 
We have found good agreement, justifying that 4 active 
layers are sufficient in the present dynamical regime. Fi- 
nally in our energetic analysis, the Ar atoms in the vicinity 
of the impact point will be emphasized. This includes an 
hemisphere around the impact point (4x4 in first layer, 
3x3 in second, and 2x2 in third). We will see that this im- 
pact hemisphere carries almost the whole energy transfer 
in the deposition process (cf. right panel of Fig. [3]). 

3 The effect of the surface 

In a first run, we investigate the influence of the mod- 
eling of the Ar substrate. The Nag center-of-mass starts 
from 15 oq above the surface with an initial velocity in 
the z direction, corresponding to a kinetic energy E^^^ of 
0.06 Ry. Experimentally, this setup should lead to stick- 
ing, if one extrapolates available experimental data with 
proper scaling laws (see the discussion below) |37| . In or- 
der to analyze the impact of the various ingredients of 
the model; Fig. [1] compares the deposition dynamics for 
three different levels of treatment, full model, dynamical 
dipoles switched off, and frozen atomic positions (while 
maintaining dynamical dipoles). The cluster is initially in 
"bottom" conflguration where the pentagon is closer to 
the surface and the top ion faces away. In all three cases, 
the first 500 fs of the dynamics proceed similar showing the 
cluster steadily approaching the surface. Then the collision 
and the subsequent evolution proceed very differently. For 
the fixed Ar atoms (lower right panel), the surface reacts 
rigidly and one observes an immediate reflection of Nag 
combined with strong internal excitations as oscillations 
of the top ion through the pentagon. For mobile Ar atoms 
without dipolc dynamics (lower left panel), there is again 
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Fig. 1. z coordinates of the Nae cluster (full curves) deposited 
on Ar384 (dotted) for an impact energy of 0.06 Ry. Three dif- 
ferent levels of treatment are shown. Upper panel: full model. 
Lower left panel: dynamical dipoles of Ar atoms switched off. 
Lower right panel: totally frozen Ar atoms. 



reflection, but now with a visible amount of excitation 
energy going to the substrate while less is left for inter- 
nal motion of the cluster. For the case of fully active Ar 
atoms (upper panel), the substrate absorbs most of the 
cluster's kinetic energy and the cluster remains tied to the 
surface, with its distance to the surface performing large 
oscillations in the attractive potential well created by the 
Ar dipoles. We thus observe a relatively soft landing of 
the Nae in "bottom" configuration. Soft means here that 
the surface and the cluster are still moderately perturbed. 
The overall comparison of the three calculations demon- 
strates the importance of a full dynamical treatment of the 
Ar surface, going beyond a mere Molecular Dynamics of 
Ar positions but also accounting for their polarizabilitics 
through time-dependent dipoles. 

The perturbation produced by the impact of the clus- 
ter on the first layer propagates straightforwardly through 
the substrate as a sound wave with about the speed of 
sound inside Ar bulk, namely around 20-30 ao/ps. When 
this sound wave reaches the 5th (and frozen) layer, it is 
reflected upwards. When coming back to the surface it 
transfers some residual momentum to the cluster. However 
that momentum transfer remains moderate and leaves the 
cluster in a captured state. That back-flow of momentum 
would not show up in that pronounced manner for an in- 
finitely deep material. We have checked that in this energy 
range and in the case of the deposit of a single Na atom, 
adding 2 more Ar layers. This indeed does not change the 
scenario qualitatively as the 2 added active Ar shells do 
absorb a rather limited amount of kinetic energy. It might 
nevertheless modify in fine the deposition threshold quan- 
titatively. However the model with frozen bottom layers 
has a realistic touch. It provides a zeroth-order simulation 
of surfaces from hard materials covered by a thin layer of 
rare gas [38]. Mind nevertheless that the case may differ 



at a quantitative level. In the case of Ar coated metals, for 
example, the underlying metal provides an extra attrac- 
tion on the deposited cluster. This may thus modify the 
details of the deposition scenario, especially at the side of 
energetic considerations and deposition thresholds. 

Altogether, Fig. [1] demonstrates the crucial role played 
by the elasticity of the surface. Ar is extremely soft and 
serves as a true stopper material for gentle deposition. The 
proper dynamical treatment of the surface also appears as 
essential. 



4 Energetic analysis 

In this section, we study the dependence on the initial ki- 
netic energy E'^^^ of Nag. This amounts to change its initial 
velocity along the z direction. Fig. [2] summarizes results 
for the time evolution viewed through spatial coordinates 
and kinetic energies. A very broad range of initial kinetic 
energies -Ej^j^, from 0.6 Ry down to 0.006 Ry, is consid- 
ered. The detailed atomic and ionic z coordinates (left 
column) show at first glance an overall similarity of all the 
different cases, except for the highest energy. There is a 
very fast stopping and quick capture of the Na cluster fol- 
lowed by persistent oscillations of the cluster distance plus 
some internal oscillations. The cluster momentum trans- 
ferred at impact propagates as a sound wave through the 
substrate. The perturbation of the substrate, of course, 
increases with the initial energy. The case of the highest 
energy differs. One sees reflection of the cluster, however 
at the price of severe destruction of the surface. It corrob- 
orates the view that Ar is an extremely efficient stopper 
material. At this point, a comparison with experiments 
performed on comparable systems |37) can be made. The 
measurements of [37] were, in fact, carried out with the 
much heavier Ag material but once deposition energies are 
properly scaled, our results are in qualitative agreement. 
Note that sufficiently high deposition energies were not at- 
tained in the experiments of [37j to access the destructive 
regime. In our simulations, the threshold for destruction 
of the substrate seems to lie between 0.05 and 0.1 Ry per 
metal atom. It should finally be noted that for all cases 
of non-destructive deposition, the course of the process 
does not depend much on the initial projectile velocity. 
The same result has been found in studies of deposition 
on large, but finite, Ar clusters [39] . 

The kinetic energies arc shown in the second column of 
Fig. [21 The kinetic energy of the Na cluster increases be- 
fore contact due to the long-range polarization interaction 
which is attractive. The additional acceleration depends 
on the initial velocity. Slower velocities allow for more en- 
ergy gain, in relative value, since the cluster moves for a 
longer time in the attractive regime. The kinetic energy is 
multiplied almost by a factor 8 for the lowest E^^^ = 0.006 
Ry whereas only insignificant acceleration occurs for the 
fastest collision. At time of impact, there emerges a very 
fast and almost complete energy transfer from the Na clus- 
ter to the substrate. In less than 0.5 ps, the Ar carry 
away almost the whole Na cluster kinetic energy and only 
very little residual kinetic energy is left to the Na cluster. 
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Fig. 2. Time evolution of deposition of Nae on Ar3g4 for five different initial kinetic energies E^^^ as indicated. Left column: 
z coordinates of Na"*" ions (heavy solid line) and of Ar atoms (faint dashed). Right column: kinetic energies for the Na cluster 
(solid lines) and for the Ar system (dashed). 



Longer times are needed for the final relaxation processes. 
Note also tlie revival of the cluster kinetic energy a few ps 
after impact. This is due to the come back of the reflected 
wave in the substrate, as discussed above. 

It is instructive to analyze energy transfers in terms 
of temperature. To that end, we calculate the intrinsic ki- 
netic energy by subtracting the contribution from the 
center-of-mass motion which is particularly relevant for 
the Na cluster. The kinetic temperature is then defined as 
Tkin = 2E'^^J3N, where iV = 6 for the cluster and 384 for 
the substrate. The time evolution of the kinetic tempera- 



ture is plotted in Fig. [3l The first column shows Tkin for 
the cluster. The initial phase is purely ccntcr-of-mass mo- 
tion without intrinsic excitation. The temperature jumps 
at the moment of impact due to the large perturbation of 
all constituents in the impact zone. The jump ends close to 
the final temperature and there remains a slow and mod- 
erate relaxation to thermal equilibrium. The right column 
of Fig. [3] shows the temperatures for the substrate in two 
ways, taken over the atoms in the vicinity of the impact 
point (the "impact hemisphere" as explained at the end of 
section [2]) or over all Ar atoms outside the impact hemi- 
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Fig. 3. Time evolution of deposition of Nag on Ar384 for five different initial kinetic energies -Ej^in as indicated (Figure [2] 
continued). Left column: kinetic temperature of Na+ ions. Right column: kinetic temperature of the Ar substrate taken over 
the impact hemisphere (dashed) denoted as "local" and over all remaining atoms (dotted) denoted "all". 



sphere. The differences between these both temperatures 
are huge. The initial energy transfer goes preferably to the 
region around the impact point. The kinetic temperature 
of the impact hemisphere shows recurrent bumps. They 
are related to slow oscillations within the substrate (see 
first column of Fig. [2]) and associated energy exchange be- 
tween potential and kinetic energy. The average tempera- 
ture relaxes slowly to that of the total system. The ampli- 
tude of the oscillations and the relaxation time strongly 
depend on the initial energy, ranging from about 5 ps for 
the weakly excited cases to outside our simulation time 
for the heftier processes. Note that the average tempera- 
ture for the highest energy is above the melting point of 
about 84 K [iOlUT] . indicating once more the destruction 



of the substrate. The electronic excitation during the colli- 
sion amounts to small dipole oscillation with an amplitude 
of about 0.05 ao and related energy content of about 6.8 
meV. 

In the second column of Fig. [2l we had tracked the 
kinetic energy transfer from the cluster to the surface as 
a function of time. It is also interesting to study the fi- 
nal repartition of the initial energy. For that purpose, the 
values of the ionic kinetic, the atomic potential and the 
atomic kinetic energies after 6 ps are recorded and are nor- 
malized to the maximum kinetic energy, -Ekin™^'^, reached 
before impact. The thus obtained energy ratios are plot- 
ted as function of E^^^^ in Fig. 21 The figure shows that 
the energy share left for the Na cluster increases with in- 
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creasing deposition energy. Nonetheless, the energy loss 
at the side of the Na cluster is dramatic, even for the 
most violent case. For the Ar substrate, we see an equal 
share between potential and kinetic energies, except for 
the highest initial energy. The gain in potential energy is 
related to the moderate spatial rearrangements in the Ar 
substrate. For E'^jjj = 0.6 Ry, the kinetic energy is 2/3 
times larger than the potential one. Indeed the emitted 
Ar in this case are blown away with a substantial velocity 
(see upper left panel of Fig. [2]). 
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Fig. 4. Partitioning of final energies over cluster kinetic en- 
ergy, substrate kinetic, and substrate potential energy normal- 
ized to initial energy. 



5 Geometrical effects 

There are different possible initial geometries (orientation, 
position relative to surface) which may produce differ- 
ent reaction pattern (for an example, see the case of Nag 
on NaCl [TFIIT^ V All results reported above are obtained 
from the starting configuration displayed on the left col- 
umn of Fig. [5l denoted by "bottom" : the top ion faces 
away from the surface and the cluster center axis is placed 
above an interstitial position of the first Ar layer. Two 
other configurations have been studied. The "centered" 
one (middle column of Fig. [5]) is similar to the "bottom" 
configuration but with the axis exactly above an Ar atom 
of the first layer. The "on top" configuration (right column 
of Fig. [5]) is obtained from the "bottom" configuration by 
reversing the top ion to face towards the surface such that 
the top ion hits the Ar surface first. In all cases, the clus- 
ter centcr-of-mass is initially positioned 15 og above the 
first layer. 

The results for the time evolution of the z coordinates 
and the kinetic energies are shown in Fig.O At first glance, 
all three results look very similar. Interesting differences 
become apparent in the details. The average kinetic energy 
of the cluster is largest for the "bottom" configuration 
and smallest for "on top". This relates to the amount of 



Ar core repulsion presented to the cluster. For "bottom" , 
the axis lies on an interstitial position and the Na ions 
stay closer to the sites of the Ar atoms. For "centered", 
the most repulsive site points in between the Na ions of 
the pentagon, the interaction being that way weakened. 
The situation is somehow similar for "on top" where the 
now closest top ion dives into the interstitial site with 
maximum distance to the atoms. This also explains why 
the top ion does not go through the pentagon in that 
case. There is also a slight difference in the final distance 
to the surface. The "on top" configuration seems to stay 
closer which may be, again, due to a minimization of core 
repulsion in this configuration. 



6 Size and charge effects 

In the previous sections, we finally observed in most cases 
the soft landing of the Nae on the Ar surface. One may 
suspect that this is favored by the oblate shape of this clus- 
ter which provides a large contact area. It is interesting 
to test the deposition scenario for different cluster geome- 
tries. A good alternative is Nag. As is well known the 8 
valence electrons of Nag form a closed shell [42j[43] (magic 
number of electrons), which leads to a close to spherical 
(up to the impact of ions) electronic shape. The overall 
shape of Nag is, in turn, itself close to spherical and is 
thus significantly different from the especially oblate one 
of Nag. The left part of Fig. [6] shows the time evolution 
for Nag impinging on Ar3g4 surface. The scenario for Nag 
is very similar to what we had observed for Nag, with effi- 
cient stopping, sudden energy transfer, and final capture. 
One even observes less excitation and deformation of the 
Nag cluster than in the Nag cases. This effect is probably 
due to the fact that Nag is electronically magic and thus 
especially stable. 

The right part of Fig.[6]shows the result for the charged 
cluster Na^. It has also the magic electron number 8 and 
is near to spherical shape. While the overall bouncing sce- 
nario is qualitatively similar (with sound wave bounce), 
the net result is quantitatively different. Differences show 
up already in the initial phase. The finite net charge yields 
stronger attractive polarization interaction with the sub- 
strate and thus acquires much more acceleration. The thus 
higher impact energy and the stronger interaction produce 
a more violent collision. Most remarkable is the closer at- 
tachment to the surface which, again, is due to the much 
stronger polarization interaction. The bottom of the Na^ 
cluster is even merging into the surface layer. This small 
example indicates that the deposition of charged clusters 
is a most interesting topic which deserves further studies, 
also in direct relation with experiments [571IH] . 



7 Conclusions 

In this paper, we have investigated the deposition dynam- 
ics of small Na clusters on an Ar(OOl) surface. To that end, 
we have employed a hierarchical model treating the Na 
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Fig. 5. Time evolution of z coordinates (lower panels) and kinetic energies (middle panels) for three different collision geometries 
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Fig. 6. Time evolution of deposition for Nag (left part) and Natj" on Ar3g4 surface at impact energy of 0.001 Ry per Na ion. 
Shown are the z coordinates (lower) as well as the kinetic energies of Na ions and Ar atoms (upper). 
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cluster in full detail by time-dependent density-functional 
theory, the Na"*" ions by classical dynamics, and the Ar 
atoms also at a classical level with position aird dipole 
polarization as dynamical variables. We have studied in 
detail the influence of initial conditions on the deposition 
process, varying substrate activity, impact energy, coUi- 
sional geometry, and size as well as charge of the Na clus- 
ter. 

We have found that the gross features of deposition 
dynamics are much the same in all cases. The Ar system 
acts always as a very efficient and smooth stopper which 
absorbs almost all (more than 90%) of the impact energy 
and thus yields capture over a broad range of energies. 
Cluster reflection, which is typical for hard surfaces, can 
be forced only at such high impact energies that the sur- 
face is destroyed. The impinging Na cluster is stopped at 
first impact and its energy is transferred to the Ar sys- 
tem in less than 0.5 ps. The transfer is almost complete 
leaving less than 10% of the energy to the cluster. The 
transferred energy is distributed also very quickly over all 
Ar atoms, propagating like a sound wave with speed of 
sound through the medium. The full thermalization over 
all ionic and atomic degrees of freedom proceeds on a slow 
scale with time of 5 ps and more (increasing with increas- 
ing impact energy). The coUisional geometry (initial ori- 
entation and horizontal position relative to the surface) 
has no effect on the overall dynamical pattern, but deter- 
mines subtle details of the captured state as, e.g., remain- 
ing cluster energy or distance to the surface. Cluster size 
and shape seem to make little difference. Cluster charge 
however can change a lot. It enhances overall excitation as 
well as binding of the cluster to the surface and produces 
a much tighter final configuration. 
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A The Na-Ar energy functional in detail 

The degrees of freedom of the model are the wavefunc- 
tions of valence electrons of the metal cluster, {(y9„(r), n — 
l...Nc\}, the coordinates of the cluster's Na+ ion cores, 
{R/, / = l...Afion}, of the Ar atoms cores Ar'5+), {Ra, a = 
l...AfAi }, and of the Ar valence clouds, {R'a, a = 1...7VAr}- 
^From the given total energy, the corresponding equations 
of motion are derived in a standard manner by variation. 
This leads to the (time-dependent) Kohn-Sham equations 
for the one single-particle wavefunctions (pn (r) of the clus- 
ter electrons, and Hamiltonian equations of motion for the 
other three degrees of freedom, thus treated by classical 
molecular dynamics (MD). For the valence cluster elec- 
trons, we use a density functional theory at the level of the 
time-dependent local-density approximation (TDLDA) , aug 
mented with an average-density self-interaction correction 



(ADSIC) [45]. The density of these electrons is given nat- 
urally as defined in mean-field theories and reads /Ooi(r) = 
l'/'n('')l • atom is described by two constituents 

with opposite charge, positive Ar core and negative Ar 
valence cloud, which allows a correct description of polar- 
ization dynamics. In order to avoid singularities, we as- 
sociate a smooth (Gaussian) charge distribution to both 
constituents having width CAr of the order of the p shell 
"size" in Ar atoms, in the spirit of [33] : 



PAr,a(r) 



eQ 



exp 



(r-RJ 



exp 



(r-R:,)^ 



'Ar 



(1) 



The corresponding Coulomb potential exerted by the Ar 
atoms is related to the charge distribution ([1]) by the Pois- 
son equation, and reads: 



crf(|r-RJ/aA..) 
|r-Ral 
CTfi\r-K\/a^^) 



r-R' 



(2) 



where erf(r') = dec e ^ stands for the error func- 

tion. As for the Na+ ions, their dynamical polarizability 
is neglected and we treat them simply as charged point 
particles. 

The total energy of the system is composed as: 



-Etotal = E] 



Naclustcr 



(3) 



The energy of the Na cluster E'Naciustcr consists out of 
TDLDA (with SIC) for the electrons, MD for ions, and 
a coupling of both by soft, local pseudo-potentials, for 
details see [?T1[501[^ . The Ar system and its coupling to 
the cluster are described by 



a 



a<a' 



drpAr,a(r)FlP°'V) 



-J r(corc) ^-p -p 



I,a 

drp„(r)E K°i^« + ^ci.Ar(|r-Ra|) 

a 

''y- r(/drfa(r)pci(r))^ 



(4) 



(5) 



VdW — 



fa(r) = V 



2 ^ "L 

a 

^ j drf,(r)2p,i(r)' 
erfdr-RJ/a^J 



Ir-R„ 



(6) 
(7) 



The Van der Waals interaction between cluster elec- 
trons and Ar dipoles is written in Eq. ([6]) as a correlation 
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from the dipolc excitation in the Ar atom coupled with a 
dipole excitation in the cluster. We exploit that the plas- 
mon frequency wmic is far below the excitations in the Ar 
atom. This simplifies the term to the variance of the dipole 
operator in the cluster, using the regularized dipole opera- 
tor fa, defined in Eq.([7]), corresponding to the smoothened 
Ar charge distributions p2ll46] . The full dipole variance is 
simplified in terms of the local variance. 

The interaction of one Ar atom with the other con- 
stituents (Ar atoms, Na"'' ions, cluster electrons) results 
from the balance between a strong repulsive core poten- 
tial that falls off exponentially and an equally strong at- 
traction from dipole polarizability. The (most important) 

polarization potential, V^°^\ is described by a valence 
electron cloud oscillating against the Ar core ion. Its pa- 
rameters are the effective charge of valence cloud Q, the 
effective mass of valence cloud ttiav — QtitIci, the restor- 
ing force for dipoles fcAr, and the width of the core and 
valence clouds UAr- The Q and fcAr are adjusted to repro- 
duce experimental data on dynamical polarizability a^i (w) 
of the Ar atom at low frequencies, namely the static limit 
a£)(w = 0) and the second derivative of a'^{J' =0). The 
width aAr is determined consistently such that the restor- 
ing force from the folded Coulomb force (for small dis- 
placements) reproduces the spring constant /sAr- 

The short range repulsion is provided by the various 
core potentials. For the Ar-Ar core interaction in Eq. (|4]), 
we employ a Lennard- Jones type potential with parame- 
ters reproducing binding properties of bulk Ar : 



^Irr^^) = e'^Ar (i?Ar/i?)'' - (i?Ar/i?)' 



(8) 



The Na-Ar core potential V{^^ in Eq. ^ is chosen ac- 
cording to [34], within properly avoiding double counting 
of the dipole polarization-potential, hence the following 
form : 



^Na,Ar(-R) 



A 



Na- 



-/3Na-R 

R 



Na,6 



2 CkAr „ „ erf(i?/\/2crAr) 

e ±1, • Vr 

2i?3 R 



Na,8 
"R8~ 



(9) 



Finally the pseudo-potential Woi.Ar in Eq- © for the electron- 
Ar core repulsion has been modeled according to the pro- 
posal of [33] : 



Vt^ci,Ar(r) = 



(10) 



The various contributions are calibrated from indepen- 
dent sources, with a final fine tuning to the NaAr dimer 
(bond length, binding energy, and optical excitation spec- 
trum) modifying only the term WoiAr- The parameters are 
summarized in table [T] The third column of the table in- 
dicates the source for the parameters. 
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■t/(pol) 


<JAr = ^5 , KAr = -^ff- , rUAr = QArmel 

n-^ ( r, . ^ 

CTRG (^QAr 3(2^)3/2 ) 


aAr=11.08ag 
ujo = 1.755 Ry 


WolAr 


A,i=0A7 , /?ei = 1.6941 /ao , =2.2 ao 


fit to NaAr dimer |46II47| 


■[/{core) 
ArAr 


AAr=1.367 * 10"^ Ry , i?Ar=6.501 ao 


fit to bulk Ar 


''^ArNa 


/3Na= 1.7624 ajY^ , QNa= 1.815 ao , ^Na= 334.85 
CNa,6= 52.5 ao , CNa,8= 1383 ao 


after [34] 



Table 1. Parameters for the various model potentials. 
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